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ABSTRACT 

We present a description of the pipeline used to calibrate the Planck Low Rrequency Instrument (LRI) timelines into thermodynamic temperatures 
for the Planck 2015 data release, covering four years of uninterrupted operations. As in the 2013 data release, our calibrator is provided by the spin- 
synchronous modulation of the cosmic microwave background dipole, but we now use the orbital component, rather than adopting the Wilkinson 
Microwave Anisotropy Probe (WMAP) solar dipole. This allows our 2015 LRI analysis to provide an independent Solar dipole estimate, which 
is in excellent agreement with that of HRI and within Icr (0.3 % in amplitude) of the WMAP value. This 0.3 % shift in the peak-to-peak dipole 
temperature from WMAP and a global overhaul of the iterative calibration code increases the overall level of the LRI maps by 0.45 % (30 GHz), 
0.64 % (44 GHz), and 0.82 % (70 GHz) in temperature with respect to the 2013 Planck data release, thus reducing the discrepancy with the power 
spectrum measured by WMAP. We estimate that the LRI calibration uncertainty is now at the level of 0.20 % for the 70 GHz map, 0.26 % for the 
44 GHz map, and 0.35 % for the 30 GHz map. We provide a detailed description of the impact of all the changes implemented in the calibration 
since the previous data release. 

Key words, cosmic microwave background - instrumentation: polarimeters - methods: data analysis 


1. Introduction 

This paper, one of a set associated with the 2015 release of 
data from the Planck^ mission, describes the techniques we em¬ 
ployed to calibrate the voltages measured by the Low Frequency 


* Corresponding author: Maurizio Tomasi, mailto imaurizio. 
tomasi@uniini. it. 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 


Instrument (LFI) radiometers into a set of thermodynamic tem¬ 
peratures, which we refer to as “photometric calibration.” We ex¬ 
pand here the work described in Planck Collaboration V (2014), 
henceforth Call3; we try to follow as closely as possible the 
structure of the earlier paper to help the reader understand what 
has changed between the 2013 and the 2015 Planck data re¬ 
leases. 

The calibration of both Planck instruments (for HFI, see 
Planck Collaboration VII 2015) is now based on the small 
(270yuK) dipole signal induced by the annual motion of the satel- 

and funded by Denmark, and additional contributions from NASA 
(USA). 
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lite around the Sun - the orbital dipole, which we derive from our 
knowledge of the orbital parameters of the spacecraft. The cali¬ 
bration is thus absolute, and not dependent on external measure¬ 
ments of the larger solar (3.35 mK) dipole, as was the case for 
Call3. Absolute calibration allows us both to improve the cur¬ 
rent measurement of the solar dipole (see Sect. 5), and to transfer 
Planck's calibration to various ground-based instruments (see, 
e.g., Perley & Butler 2015) and other cosmic microwave back¬ 
ground (CMB) experiments (e.g., Louis et al. 2014). 

Accurate calibration of the LFI is crucial to ensure reliable 
cosmological and astrophysical results from the Planck mission. 
Internally consistent photometric calibration of the nine Planck 
frequency channels is essential for component separation, where 
we disentangle the CMB from the varoius Galactic and extra- 
galactic foreground emission processes (Planck Collaboration 
IX 2015; Planck Collaboration X 2015). In addition, the LFI 
calibration directly impacts the Planck polarization likelihood 
at low multipoles, based on the LFI 70 GHz channel, which is 
extensively employed in the cosmological analysis of this 2015 
release. Furthermore, a solid absolute calibration is needed to 
compare and combine Planck data with results from other ex¬ 
periments, most notably with WMAP. Detailed comparisons be¬ 
tween calibrated data from single LFI radiometers, between the 
three LFI frequency channels, and between LFI and HFI, allow 
us to test the internal consistency and accuracy of our calibra¬ 
tion. 

In this paper, we quantify both the absolute and relative ac¬ 
curacy in the calibration of the LFI instrument, and find an over¬ 
all uncertainty of 0.35 % (30 GHz map), 0.26% (44 GHz), and 
0.20 % (70 GHz). The level of the power spectrum near the first 
peak is now remarkably consistent with WMAP’s. Other paper 
in this Planck data release deal with the quality of the LFI cali¬ 
bration, in particular: 


- Planck Collaboration X (2015) quantifies the consistency be¬ 
tween the calibration of the LFI/HFI/WMAP channels in the 
context of foreground component separation, finding that the 
measured discrepancies among channels are of the order of 
a few tenths of a percent; 

- Planck Collaboration XI (2015) analyses the consistency 
between the LFI 70 GHz low-^ polarization map and the 
WMAP map in pixel space, finding no hints of inconsisten¬ 
cies; 

- Planck Collaboration XIII (2015) compares the estimate for 
the T and Zre cosmological parameters (reionization optical 
depth and redshift) using either LFI 70 GHz polarization 
maps or WMAP maps, finding statistically consistent values. 


To achieve calibration accuracy at the few-per-thousand level 
requires careful attention to instrumental systematic effects and 
foreground contamination of the orbital dipole. Much of this pa¬ 
per is devoted to a discussion of such effects and the means to 
mitigate them. 

In this paper we do not explicitly discuss polarization-related 
issues. Although polarization analysis is one of the most im¬ 
portant results of this data release, the calibration of the LFI 
radiometers is inherently based on temperature signals (Leahy 
et al. 2010). Estimates of the sensitivity in polarization, as well 
as the impact of calibration-related systematics on it, are pro¬ 
vided by Planck Collaboration III (2015). 



Fig. 1. Schematic of an LFI radiometer, taken from Call 3. The 
two linearized polarization components are separated by an or¬ 
thomode transducer (OMT), and each of them enters a twin ra¬ 
diometer, only one of which is shown in the figure. A first ampli¬ 
fication stage is provided in the cold (20 K) focal plane, where 
the signal is combined with a reference signal originating in a 
thermally stable 4.5-K thermal load. The radio frequency sig¬ 
nal is then propagated through a set of composite waveguides 
to the warm (300 K) back-end, where it is further amplified and 
filtered, and finally converted into a sequence of digitized num¬ 
bers by an analogue-to-digital converter. The numbers are then 
compressed into packets and sent to Earth. 

1.1. Basis of the calibration 

A schematic of the LFI pseudo-correlation receiver is shown in 
Fig. 1. We model the output voltage V(t) of each radiometer as 

V(t) = G(t) [B * (Tsky + D)](t) + M + A, (1) 

where G is the “gain” (measured in VK“^), B is the beam re¬ 
sponse, D is the thermodynamic temperature of the total CMB 
dipole signal (i.e., a combination of the Solar and orbital compo¬ 
nents, including the quadrupolar relativistic corrections), which 
we use as a calibrator, and Tsky = Tcmb + T'cai + T'other is the over¬ 
all temperature of the sky (CMB anisotropies, diffuse Galactic 
emission and other^ foregrounds, respectively) apart from D. 
Finally, M is a constant offset and A is a noise term. Note that 
in the following sections we will use Eq. (1) many times; when¬ 
ever the presence of the A term will not be important, it will be 
silently dropped. The * operator represents a convolution over 
the 471 sphere. We base our calibration on the a priori knowl¬ 
edge of the spacecraft velocity around the Sun, producing the 
orbital dipole and use the orbital dipole to accurately measure 
the dominant solar dipole component. The purpose of this pa¬ 
per is to explain how we implemented and validated the pipeline 
that estimates the “calibration constant” K = G~^ (which is used 
to convert the voltage V into a thermodynamic temperature), to 
quantify the quality of our estimate for K, and to quantify the im¬ 
pact of possible systematic calibration errors on the Planck/LFl 
data products. 

^ Within this term we include extragalactic foregrounds and all point 
sources. 
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1.2. Structure of this paper 

Several improvements were introduced in the LFI pipeline for 
calibration relative to Call3. In Sect. 2 we recall some terminol¬ 
ogy and basic ideas presented in Cal 13 to discuss the normal¬ 
ization of the calibration, i.e., what factors influence the average 
value of G in Eq. (1). Section 3 provides an overview of the new 
LFI calibration pipeline and underlines the differences with the 
pipeline described in Call3. One of the most important improve¬ 
ments in the 2015 calibration pipeline is the implementation of 
a new iterative algorithm to calibrate the data, DaCapo. Its prin¬ 
ciples are presented separately in a dedicated section. Sect. 4. 
This code has also been used to characterize the orbital dipole. 
The details of this latter analysis are provided in Sect. 5, where 
we present a new characterization of the Solar dipole. These two 
steps are crucial for the calibration of LFI. Section 6 describes a 
number of validation tests we have run on the calibration, as well 
as the results of a quality assessment. This section is divided into 
several parts: in Sect. 6.1 we compare the overall level of the cal¬ 
ibration in the 2015 LFI maps with those in the previous data re¬ 
lease; in Sect. 6.2 we provide a brief account of the simulations 
described in Planck Collaboration III (2015), which assess the 
calibration error due to the white noise and approximations in 
the calibration algorithm itself; in Sect. 6.3 we describe how un¬ 
certainties in the shape of the beams might affect the calibration; 
Sects. 6.4 and 6.5 measure the agreement between radiometers 
and groups of radiometers in the estimation of the TT power 
spectrum; and Sect. 6.6 provides a reference to the discussion of 
null tests provided in Planck Collaboration III (2015). Finally, in 
Sect. 7, we derive an independent estimate of the LFI calibration 
from our measurements of Jupiter and discuss its consistency 
with our nominal dipole calibration. 


2. Handling beam efficiency 

In this section we develop a mathematical model to relate the 
absolute level of the calibration (i.e., the average level of the raw 
power spectrum Q for an LFI map) to a number of instrumental 
parameters related to the beams and the scanning strategy. 

The beam response B{6, ip) is a dimensionless function de¬ 
fined over the An sphere. In Eq. (1), B appears in the convolution 


B * (Tsky + D) = 


+ D)(e, ip) dQ 
J,^B{0,<p)dn 


( 2 ) 


whose value changes with time because of the change of orien¬ 
tation of the spacecraft. Since no time-dependent optical effects 
are evident from the data taken from October 2009 to February 
2013 (Planck Collaboration IV 2015), we assume there is no in¬ 
trinsic change in the shape of B during the surveys. 

In the previous data release, we approximated 5 as a Dirac 
delta function (a “pencil beam”) when modeling the dipole sig¬ 
nal seen by the LFI radiometers. The same assumption has been 
used for all the WMAP data releases (see, e.g., Hinshaw et al. 
2009), as well as in the HFI pipeline (Planck Collaboration VIII 
2015). However, the real shape of B deviates from the ideal case 
of a pencil beam because of two factors: (1) the main beam is 
more like a Gaussian with an elliptical section, whose FWHM 
(full width half maximum) ranges between 13' and 33' in the 
case of the LFI radiometers; and (2) farther than 5° from the 
beam axis, the presence of far sidelobes further dilutes the sig¬ 
nal measured through the main beam and induces an axial asym¬ 


metry on B. Previous studies^ tackled the first point by applying 
a window function to the power spectrum computed from the 
maps in order to correct for the finite size of the main beam. 
However, the presence of far sidelobes might cause the presence 
of stripes in maps. For this 2015 data release, we use the full 
shape of B in computing the dipole signal adopted for the cali¬ 
bration. No significant variation in the level of the CMB power 
spectra with respect to the previous data release is expected, 
since we are basically subtracting power during the calibration 
process instead of reducing the level of the power spectrum by 
means of the window function. However, this new approach im¬ 
proves the internal consistency of the data, since the beam shape 
is taken into account from the very first stages of data process¬ 
ing (i.e., the signal measured by each radiometer is fitted with 
its own calibration signal B^^d * B>); see Sect. 6.6. The definition 
of the beam window function has been changed accordingly; see 
Planck Collaboration IV (2015). 

In Cal 13 we introduced the two quantities (pj) and 0sky as a 
way to quantify the impact of a beam window function on the 
calibration"^ and on the mapmaking process, respectively. Here 
we briefly summarize the theory, and we introduce new equa¬ 
tions that are relevant for understanding the normalization of the 
new Planck-hFl results in this data release. 

The Solar dipole D, due to the motion of the Solar System in 
the CMB rest frame, is given by 


D(x,t) = Tcmb 


(r(t)ii-m-x) 



( 3 ) 


where Tcmb is the CMB monopole, P = i;/c is the velocity of 
the spacecraft, and y = (1 - Each radiometer measures 

the signal D convolved with the beam response B, according to 
Eq. (2); therefore, in principle, each radiometer has a different 
calibration signal. Under the assumption of a Dirac delta shape 
for B, Cal 13 shows that the estimate of the gain constant G is 
related to the true gain G by the formula 

GP“ = G(l-/s,)(l+M (4) 


where 



( 5 ) 


is the fraction of power entering the sidelobes (i.e., along direc¬ 
tions farther than 5° from the beam axis), and 


dtB^i >1^ D 

^r^main * D 


( 6 ) 


is a time-dependent quantity that depends on the shape of B = 
^main + ^si and its decomposition into a “main” (6 < 5°) and 
“sidelobe” part, on the signal D, and on the scanning strategy 


^ Apart from the use of appropriate window functions (e.g.. Page 
et al. 2003), the WMAP team implemented a number of other correc¬ 
tions to further reduce systematic errors due to the non-ideality of their 
beams. In their first data release, the WMAP team estimated the con¬ 
tribution of the Galaxy signal picked up through the sidelobes at the 
map level (Barnes et al. 2003) and then subtracted them from the maps. 
Starting from the third year release, they estimated a multiplicative cor¬ 
rection, called the “recalibration factor,” assumed constant throughout 
the survey, by means of simulations. This constant accounts for the side- 
lobe pickup and has been applied to the TODs (Jarosik et al. 2007). The 
deviation from unity of this factor ranges from 0.1 % to 1.5 %. 

^ The definition of (po provided in Call3 was not LFI-specific: it can 
be applied to any experiment that uses the dipole signal for the calibra¬ 
tion. 
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Fig. 2. Quantities used in the determination of the value of 0sky (Eq* 8) for radiometer LFI-27M (30 GHz) during a short time span 
(2 min). Panel A\ the quantities 5mb * Tqmb and ^si * T'cmb are compared. Fluctuations in the latter term are much smaller than 
those in the former. Panel B\ the quantity * Tcmb shown in the previous panel is replotted here to highlight the features in its 
tiny fluctuations. The fact that the pattern of fluctuations repeats twice depends on the scanning strategy of Planck, which observes 
the sky along the same circle many times. Panel C: Value of 0sky calculated using Eq. (8). There are several values that diverge to 
inflnity, which is due to the denominator in the equation going to zero. Panel D\ distribution of the values of 0sky plotted in panel C. 
The majority of the values fall around the number +0.02 %. 


because of the time dependence of the stray light; the notation dt 
indicates a time derivative. Once the timelines are calibrated, tra¬ 
ditional mapmaking algorithms approximate^ 5 as a Dirac delta 
(e.g., Hinshaw et al. 2003; Jarosik et al. 2007; Keihanen et al. 
2010), thus introducing a new systematic error. In this case, the 
mean temperature f sky of a pixel in the map would be related to 
the true temperature Tsky by the formula 

r3ky=fP“(l-.^sky+.^D), (7) 

which applies to timelines and should be considered valid only 
when considering details on angular scales larger than the width 
of the main beam. Cal 13 deflnes the quantity 0sky using the fol¬ 
lowing equation: 

^sky 
^sky 

See Fig. 2 for an example showing how 0sky is computed. 

In this 2015 Planck data release, we take advantage of our 
knowledge of the shape of B to compute the value of Eq. (2) and 
use this as our calibrator. Since the term B * Tsky — (^main + 5sl)* 
Tsky is unknown, we apply the following simpliflcations: 

^ Keihanen & Reinecke (2012) provide a deconvolution code that can 
be used to produce maps potentially free of this effect. 


1. we apply the point source and 80 % Galactic masks (Planck 
Collaboration ES 2015), in order not to consider the ^main * 
Tsky term in the computation of the convolution; 

2. we assume that B^\ * Tsky ~ ^si * Tq^\ and subtract it from 
the calibrated timelines, using an estimate for Tq^\ com¬ 
puted by means of models of the Galactic emission (Planck 
Collaboration IX 2015; Planck Collaboration X 2015). 

The result of such transformations is a new timeline Under 
the hypothesis of perfect knowledge of the beam B and of the 
dipole signal D, these steps are enough to estimate the true cali¬ 
bration constant without bias^ (unlike Eq. 4): 

= G, (9) 

which should be expected; since no systematic effects caused 
by the shape of B affect the estimate of the gain G. In order 
to see how Eq. (7) changes in this case, we write the measured 
temperature f sky as 

^sky = B ^ Tsky + M = ^main * ^sky + ^sl * ^sky + (10) 


^ It is easy to show this analytically. Alternatively, it is enough to note 
that considering the full 4?: beam makes /si = 0, and is identically 
zero because there are no “sidelobes” falling outside the beam. With 
these substitutions, Eq. (4) becomes Eq. (9). 


0sky — 


^si * T, 


sky 


^main * Tsky 


\ ^ B^\ * Tsky 
/ ^sky 
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Since in this 2015 data release we remove B^\ * Toai, the con¬ 
tribution of the pickup of Galactic signal through the sidelobes 
(Planck Collaboration II 2015), the equation can be rewritten as 

“ ^main * T'sky + ^sl * (T'cMB + Mother) + M. (11) 

If we neglect details at angular scales smaller than the main beam 
size, then 

^main * T'sky ~ (1 “ /si) T'sky, (12) 

SO that 


TX = (1 - /si) ^’sky + Bsl * {TcMB + T’other) + M. (13) 
We modify Eq. (8) in order to introduce a new term 
^sl * (T'cMB + T'other) 


/sky 


* sky 


solving for Tsky, Eq. (13) can be rewritten as 

1 -( 


T^sky = t: 


47T 


^sky 


l-/sl 


+ To, 


(14) 


(15) 


where Tq = M/{1 - /si) is a constant offset that is of little rel¬ 
evance for pseudo-differential instruments like LEI. Eq. (15) is 
the equivalent of Eq. (7) in the case of a calibration pipeline 
that takes into account the 47r shape of B, as is the case for the 
Planck-LFl pipeline used for the 2015 data release. 

Since one of the purposes of this paper is to provide a quan¬ 
titative comparison of the calibration of this Planck data release 
with the previous one, we provide now a few formulae that quan¬ 
tify the change in the average level of the temperature fluctua¬ 
tions and of the power spectrum between the 2013 and 2015 re¬ 
leases. The variation in temperature can be derived from Eq. (7) 
and Eq. (15): 


7-2015 

sky 

7-2013 

sky 


(1 /sky + /d)( 1 /si) 


1 -/' 


sky 


1 /si /sky + < 


1-/: 


(16) 


sky 


If we consider the ratio between the power spectra and 
the quantity becomes 


r 1 - /si - /sky + /d 1 

1 ] ■ 


(17) 


In Sect. 6.1 we will provide quantitative estimates of /si, /d, /sky, 
and as well as the ratios in Eq. (16) and in Eq. (17). 


3. The calibration pipeiine 

In this section we briefly describe the implementation of the cal¬ 
ibration pipeline. Readers interested in more detail should refer 
to Planck Collaboration II (2015). 

Evaluating the calibration constant K (see Eq. 1) requires us 
to fit the timelines of each radiometer with the expected signal 
D induced by the dipole as Planck scans the sky. This process 
provides the conversion between the voltages and the measured 
thermodynamic temperature. 

As discussed in Sect. 2, we have improved the model used 
for D, since we are now computing the convolution of D with 
each beam B over the full An sphere. Moreover, we are consid¬ 
ering the B^\ * Toai term in the fit, in order to reduce the bias 
due to the pickup of Galactic signal by the beam far sidelobes. 


The model of the dipole D now includes the correct^ quadrupo- 
lar corrections required by special relativity. The quality of the 
beam estimate B has been improved as well: we are now us¬ 
ing all the seven Jupiter transits observed in the full 4-year mis¬ 
sion, and we account for the optical effects due to the varia¬ 
tion of the beam shape across the band of the radiometers. It 
is important to underline that these new beams do not follow the 
same normalization convention as in the first data release (now 
/ B(0,^)dQ * 1), as numerical inaccuracies in the simulation 
oi the An beams cause a loss of roughly 1 % of the signal en¬ 
tering the sidelobes^: see Planck Collaboration IV (2015) for a 
discussion of this point. 

As was the case in the 2013 data release, the calibration con¬ 
stant K is estimated once per each pointing period, i.e., the pe¬ 
riod during which the spinning axis of the spacecraft holds still 
and the spacecraft rotates at a constant spinning rate of 1 /60 Hz. 
The code used to estimate K, named DaCapo, has been com¬ 
pletely rewritten; it is able to run in two modes, one of which (the 
so-called unconstrained mode) is able to produce an estimate of 
the solar dipole signal, and the other one (the constrained mode) 
which requires the solar dipole parameters as input. We have 
used the unconstrained mode to assess the characteristics of the 
solar dipole, which have then been used as an input to the con¬ 
strained mode of DaCapo for producing the actual calibration 
constants. 

We smooth the calibration constants produced by DaCapo by 
means of a running mean, where the window size has a variable 
length. That length is chosen so that every time there is a sudden 
change in the state of the instrument (e.g., because of a change in 
the thermal environment of the front-end amplifiers) that discon¬ 
tinuity is not averaged out. However, this kind of filter removes 
any variation in the calibration constants, whose timescale is 
smaller than a few weeks. One example of this latter kind of fluc¬ 
tuation is the daily variation measured in the radiometer back¬ 
end gains during the first survey, which was caused by the con¬ 
tinuous turning on-and-off of the transponder^ while sending the 
scientific data to Earth once per day. In order to keep track of 
such fluctuations, we have estimated the calibration constants K 
using the signal of the 4.5 K reference load in a manner similar 
to that described in Cal 13 under the name of “4K calibration”, 
and we have added this estimate to the DaCapo gains after hav¬ 
ing applied a high-pass Alter to them, as shown in Eig. 5. Details 
about the implementation of the smoothing Alter are provided in 
Appendix A. 

Once the smoothing Alter has been applied to the calibra¬ 
tion constants K, we multiply the voltages by K in order to con¬ 
vert them into thermodynamic temperatures, and we remove the 
term B ^ D + B^\^ Tsky from the result, thus removing the dipole 
and the Galactic signal captured by the far sidelobes from the 
data. The value for Tsky has been taken from a sum of the fore¬ 
ground signals considered in the simulations described in Planck 
Collaboration IX (2015); refer to Planck Collaboration II (2015) 
for further details. 


^ Because of a bug in the implementation of the pipeline, the previ¬ 
ous data release had a spurious factor that led to a residual quadrupolar 
signal of ~ 1.9 fjK, as described in Call3. 

^ This loss was present in the beams used for the 2013 release too, but 
in that case we applied a normalization factor to B. The reason why we 
removed this normalization is that it had the disadvantage of uniformly 
spreading the 1 % sidelobe loss over the whole An sphere. 

^ This operating mode was subsequently changed and the transpon¬ 
der has been kept on for the remainder of the mission starting from 272 
days after launch, thus removing the origin of this kind of gain fluctua¬ 
tions. 
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Fig. 3. Diagram of the pipeline used to produce the LFI frequency maps in the 2015 Planck data release. The grey ovals represent 
input/output data for the modules of the calibration pipeline, which are represented as white boxes. The product of the pipeline is a 
set of calibrated timelines that are passed as input to the mapmaker. 


4. The calibration aigorithm 

DaCapo is an implementation of the calibration algorithm we 
have used in this data release to produce an estimate of the cal¬ 
ibration constant K in Eq. (1). In this section we describe the 
model on which DaCapo is based, as well as a few details of its 
implementation. 

4.1. Unconstrained algorithm 

Let Vi be the ith sample of an uncalibrated data stream, and 
k(i) the pointing period to which the sample belongs. Following 
Eq. (1) and assuming the usual mapmaking convention of scan¬ 
ning the sky, using a pencil beam, we model the uncali¬ 
brated time stream as 


V.- = Gkii)(Ti + B * A) + bkii) + A, (18) 

where we write B ^ Di = {B ^ D)i and use the shorthand notation 
Ti = (Tsky)/. The quantity Gk is the unknown gain factor for ^th 
pointing period, nt represents white noise, and bk is an offset^^ 
which captures the correlated noise component. We denote by 
Ti the sky signal, which includes foregrounds and the CMB sky 
apart from the dipole, and hy B ^ Di the dipole signal as seen 
by the beam B. The dipole includes both the Solar and orbital 
components, and it is convolved with the full 471 beam. The beam 

The offset absorbs noise at frequencies lower than the inverse of 
the pointing period length (typically 40 min). The process of coadding 
scanning rings effciently reduces noise at higher frequencies. We treat 
the remaining noise as white. 


convolution is carried out by an external code, and the result is 
provided as input to DaCapo. 

The signal term is written with help of a pointing matrix P 
as 

Ti = y]Pipmp. (19) 

P 

Here P is a pointing matrix which picks the time-ordered sig¬ 
nal from the unknown sky map m. The current implementa¬ 
tion takes into consideration only the temperature component. 
In radiometer-based calibration, however, the polarisation signal 
is partly accounted for, since the algorithm interprets as tem¬ 
perature signal whatever combination of the Stokes parameters 
(/, 2, U) the radiometer records. In regions that are scanned in 
one polarisation direction only, this gives a consistent solution 
that does not induce any error on the gain. A small error can be 
expected to arise in those regions where the same sky pixel is 
scanned in vastly different directions of polarisation sensitivity. 
The error is proportional to the ratio of the polarisation signal 
and the total sky signal, including the dipole. 

We determine the gains by minimising the quantity 

= y 4(y. _ ymod)2^ ^20) 

^ cr 


where 


T /mod 

L- = Gk{i) 


^ Pipmp -\- B ^ Di 




’k(0- 


( 21 ) 


\ p / 

and cr^ is the white noise variance. The unknowns of the model 
are m, G, b, and n (while ee assume that the beam B is perfectly 
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Fig. 4. Variation in time of a few quantities relevant for calibration, for radiometer LFI-21M (70 GHz, left) and LFI-27M (30 GHz, 
right). Grey/white bands mark complete sky surveys. All temperatures are thermodynamic. Panel A\ calibration constant K estimated 
using the expected amplitude of the CMB dipole. Note that the uncertainty associated with the estimate changes with time, according 
to the amplitude of the dipole as seen in each ring. Panel B: expected peak-to-peak difference of the dipole signal (solar + orbital). 
The shape of the curve depends on the scanning strategy of Planck, and it is strongly correlated with the uncertainty in the gain 
constant (see panel A). Note that the deepest minima happen during Surveys 2 and 4; because of the higher uncertainties in the 
calibration (and the consequent bias in the maps), these surveys have been neglected in some of the analyses in this Planck data 
release (see e.g., Planck Collaboration XIII 2015). Panel C: the calibration constants K used to actually calibrate the data for this 
Planck data release are derived by applying a smoothing filter to the raw gains in panel A. Details regarding the smoothing filter are 
presented in Appendix A. 


known). The dipole signal D and pointing matrix P are assumed 
to be known. 

To reduce the uncertainty that arises from beam effects and 
subpixel variations in signal, we apply a galactic mask and in¬ 
clude in the sum in Eq. (20) only those samples that fall outside 
the mask. 

Since Eq. (21) is quadratic in the unknowns, the minimisa¬ 
tion of;^^ requires iteration. To linearise the model we first rear¬ 
range it as 


P 

p 

+ - m°)] + bkn). 


Here and mP are the gains and the sky map from the previous 
iteration step. We drop the quadratic term in brackets and obtain 


that 


vr" = G,(,) 


B * D, + ^ Pintri 


p 

IP"^P 


■ ^^(0 ^ Pip^p + ^k{i)- 


(23) 


Here 


fhp = {nip - nPp) 


(24) 


is a correction to the map estimate from the previous iteration 
step. Eq. (23) is linear in the unknowns m, G and b. We run an 
iterative procedure, where at every step we minimizewith the 
linearized model in Eq. (23), update the map and the gains as 
wP m and ^ g, and repeat until convergence. The 
iteration is started from = 0. Thus at the first step we 

are fitting just the dipole model and a baseline Gk(i) B^Di-\-bk, and 
we obtain the first estimate for the gains. The first map estimate 
is obtained in the second iteration step. 

DaCapo solves the gains for two radiometers of a horn at the 
same time. Two map options are available. Either the radiome- 
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Dipole fit + DaCapo 
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Fig. 5. Visual representation of the algorithms used to filter the calibration constants produced by DaCapo (top-left plot; see Sect. 4). 
The example in the figure refers to radiometer LFI-27M (30 GHz) and only shows the first part of the data (roughly three surveys). 
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ters have their own sky maps, or both see the same sky. In the 
former case the calibrations become independent. 


4.1.1. Solution of the linear system 

Minimisation of yields a large linear system. The number of 
unknowns is dominated by the number of pixels in map m. It is 
possible, however, to reformulate the problem as a much smaller 
system as follows. 

We first rewrite the model using matrix notation. We com¬ 
bine the first and last terms of Eq. (23) formally into 


nil) 




ip^p 


+ bi 


kii) 




(25) 


The vector aj contains the unknowns b and G, and the matrix F 
spreads them into a time-ordered data stream. The dipole signal 
D seen by the beam B, and a signal picked from map nP, are 
included in F. 

Eq. (21) can now be written in matrix notation as 

ymodel ^ (26) 

Gains have been transferred inside matrix P, 

Pip = (27) 

Using this notation, Eq. (20) becomes 

= (V-Pm-Fa)^C;*(V-Pm-Fa), (28) 

where is the white noise covariance. 

Eq. (28) is equivalent to the usual destriping problem of map¬ 
making (Planck Collaboration VI 2015), only the interpretation 
of the terms is slightly different. In place of the pointing matrix 
P we have P, which contains the gains from previous iteration 
step, and a contains the unknown gains beside the usual baseline 
offsets. 

We minimize Eq. (28) with respect to m, insert the result 
back into Eq. (28), and minimize with respect to a. The solution 
is identical to the destriping solution 

a = (FTc-'ZF)-'FTc-‘Zy, (29) 

where 

Z = I - P(pTc-‘ P)-‘ Fc-‘. (30) 

We use a hat to indicate that a is an estimate of the true a. We are 
here making use of the sparse structure of the pointing matrix, 
which allows us to invert matrix P through non-iterative 

methods. Eor a detailed solution of an equivalent problem in 
mapmaking, see Keihanen et al. (2010) and references therein. 
The linear system in Eq. (29) is much smaller than the original 
one. The rank of the system is of the order of the number of 
pointing periods, which is 44 070 for the full four-year mission. 

Eq. (29) can be solved by conjugate gradient iteration. The 
map correction is obtained as 


m = (FC„-‘F)-'FQ‘(V-Fa). (31) 

Matrix P^C~^P is diagonal, and inverting it is a trivial task. 

A lower limit for the gain uncertainty, based on radiometer 
white noise only, is given by the covariance matrix 

Ca = (F^C„-‘F)-'. (32) 


4.2. Constrained algorithm 

4.2.1. Role of the Solar dipole 

The dipole signal is a sum of the Solar and orbital contributions. 
The Solar dipole can be thought of as being picked from an ap¬ 
proximately^^ constant dipole map, while the orbital component 
depends on beam orientation and satellite velocity. The latter can 
be used as an independent and absolute calibration. As we will 
discuss in Sect. 5, this has allowed us to determine the ampli¬ 
tude and direction of the Solar dipole and decouple the Planck 
absolute calibration from that of WMAR 

The Solar dipole can be interpreted either as part of the 
dipole signal B ^ D ox part of the sky map m. This has impor¬ 
tant consequences. The advantage is that we can calibrate us¬ 
ing only the orbital dipole, which is better known than the solar 
component and can be measured absolutely (it only depends on 
the temperature of the CMB monopole and the velocity of the 
Planck spacecraft). When the unconstrained DaCapo algorithm 
is run with erroneous dipole parameters, the difference between 
the input dipole and the true dipole simply leaks into the sky 
map m. The map can then be analysed to yield an estimate for 
the Solar dipole parameters. 

The drawback from the degeneracy is that the overall gain 
level is weakly constrained, since it is determined from the or¬ 
bital dipole alone. In the absence of the orbital component, a 
constant scaling factor applied to the gains would be fully com¬ 
pensated by an inverse scaling applied to the signal. It would 
then be impossible to determine the overall scaling of the gain. 
The orbital dipole breaks the degeneracy, but leaves the overall 
gain level weakly constrained compared with the relative gain 
fluctuations. 

The degeneracy is not perfect, since the signal seen by a ra¬ 
diometer is modified by the beam response B. In particular, a 
beam sidelobe produces a strongly orientation-dependent signal. 
This is however, a small correction to the full dipole signal. 


4.2.2. Dipole constraint 

Because of the degeneracy between the overall gain level and the 
map dipole, it makes sense to constrain the map dipole to zero. 
Eor this to work, two conditions must be fulfilled: 1) the Solar 
dipole must be known, and 2) the contribution of foregrounds 
(outside the mask) to the dipole of the sky must either be negli¬ 
gible, or it must be known and included in the dipole model. 

In the following we assume that both the orbital and the Solar 
dipole are known. We aim at deriving a modified version of the 
DaCapo algorithm, where we impose the additional constraint 
mjm = 0. Here m/) is a a map representing the Solar dipole 
component. We are thus requiring that the dipole in the direction 
of the Solar dipole is completely included in the dipole model D, 
and nothing is left for the sky map. Note that iri/) only includes 
the pixels outside the mask. 

It turns out that condition m^m = 0 alone is not sufficient, 
since there is another degeneracy in the model that must be taken 
into account. The monopole of the sky map is not constrained 
by data, since it cannot be distinguished from a global noise 
offset bk=constSint. It is therefore possible to satisfy the condi¬ 
tion mjm = 0 by adjusting simultaneously the baselines and the 

It is not exactly constant, as the dipole signal is B * D. Since the 
orientation of B changes with time, any deviation from axial symmetry 
in B (ellipticity, far sidelobes...) falsifies this assumption. However, 
when convolving a large-scale signal such as the CMB dipole with the 
LFI beams, such asymmetries are a second-order effect. 
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monopole of the map, with no cost in;^^. To avoid this pitfall, 
we constrain simultaneously the dipole and the monopole of the 
map. We require m^m = 0 and l^m = 0, and combine them 
into one constraint 

m^m = 0, (33) 

where rric now is a two-column object. 

We add to Eq. (28) an additional prior term 

y2 = (V - Pm - Fa)^C:‘(V - M - Pa) 

-r m nicCj m^m 

and aim at taking to infinity. This will drive m^m to zero. 
Minimisation of Eq. (34) yields the solution 

a = {F^C-^ZF)-^F^C-^ZW, (35) 

with 

Z = I- P(M -r (36) 

where we have written for brevity 

M = P^C~^P. (37) 

This differs from the original solution (Equations 29-30) by the 
term in the definition of Z. 

Eq. (36) is unpractical due to the large size of the matrix to be 
inverted. To proceed, we apply the Sherman-Morrison formula 
and let Cz) ^ 0, yielding 

(38) 

The middle matrix m^M ^nVc is a 2x2 block diagonal matrix, 
and is easy to invert. 

Equations (35)-(38) are the basis of the constrained DaCapo 
algorithm. The system is solved using a conjugate-gradient 
method, similarly to the unconstrained algorithm. 

The map correction becomes 

m = (M + m,Co'mf)-'p2'Q'(V - Pa). (39) 

One readily sees that in fulfills the condition expressed by 
Eq. (33), and thus so does the full map m. 

The constraint breaks the degeneracy between the gain and 
the signal, but also makes the gains again dependent on the Solar 
dipole, which must be known beforehand. 

4.3. Use of unconstrained and constrained algorithms 

We have used the unconstrained and constrained versions of the 
algorithm together to obtain a self-consistent calibration, and to 
obtain an independent estimate for the Solar dipole. 

We have first run the unconstrained algorithm, using the 
known orbital dipole and an initial guess for the Solar dipole. 
The results depend on the Solar dipole only through the beam 
correction. The difference between the input dipole and the true 
dipole are absorbed in the sky map. 

We have estimated the Solar dipole from these maps 
(Sect. 5); since the Solar dipole is same for all radiometers, we 
have combined data from all the 70 GHz radiometers to reduce 
the error bars. (Simply running the unconstrained algorithm, fix¬ 
ing the dipole, and running the constrained version with same 
combination of radiometers would just have yielded the same 
solution.) 

Once we have produced an estimate of the Solar dipole, we 
reran DaCapo in constrained mode to determine the calibration 
coefficients K more accurately. 


5. Characterization of the Orbitai and Soiar dipoies 

In this section, we explain in detail how the Solar dipole was ob¬ 
tained for use in the final DaCapo run mentioned above. We also 
compare the LEI measurements with the Planck nominal dipole 
parameters, and with the WMAP values given by Hinshaw et al. 
(2009). 

5.1. Analysis 

When running DaCapo in “constrained” mode to compute the 
calibration constants K (Eq. 1), the code needs an estimate of 
the Solar dipole in order to calibrate the data measured by the 
radiometers (see Sect. 3 and especially Eig. 3), since the sig¬ 
nal produced by the orbital dipole is ten times weaker. We have 
used DaCapo to produce this estimate from the signal produced 
by the orbital dipole. We limited our analysis to the 70 GHz ra¬ 
diometric data, since this is the cleanest frequency in terms of 
foregrounds. The pipeline was provided by a self-contained ver¬ 
sion of the DaCapo program, run in unconstrained mode (see 
Sect. 4.2), in order to make the orbital dipole the only source of 
calibration, while the Solar dipole is left in the residual sky map. 

We bin the uncalibrated differenced time-ordered data into 
separated rings, with one ring per pointing period. These data 
are then binned according to the direction and orientation of the 
beam, using aHealpix^^ (Gorski et al. 2005) map of resolution 
A^side = 1024 and 256 discrete bins for the orientation angle if/. 
The far sidelobes would prevent a clean dipole from being recon¬ 
structed in the sky model map, since the signal in the timeline is 
convolved with the beam B over the full sphere. To avoid this, an 
estimate for the pure dipole is obtained by subtracting the contri¬ 
bution due to far sidelobes using an initial estimate of the Solar 
dipole, which in this case was the WMAP dipole (Hinshaw et al. 
2009). We also subtract the orbital dipole at this stage. The bias 
introduced by using a different dipole is of second order and is 
discussed in the section on error estimates. 

DaCapo builds a model sky brightness distribution that is 
used to clean out the polarized component of the CMB and fore¬ 
ground signals to leave just noise, the orbital dipole, and far- 
sidelobe pickup. This sky map is assumed to be unpolarized, 
but since radiometers respond to a single linear polarization, the 
data will contain a polarized component, which is not compati¬ 
ble with the sky model and thus leads to a bias in the calibration. 
An estimate of the polarized signal, mostly CMB £’-modes and 
some synchrotron, needs to be subtracted from the timelines. To 
bootstrap the process we need an intial gain estimate, which is 
provided by DaCapo constrained to use the WMAP dipole. We 
then used the inverse of these gains to convert calibrated po¬ 
larization maps from the previous LEI data release (which also 
used WMAP dipole calibration) into voltages and unwrap the 
map data into the timelines using the pointing information for 
position and Foresight rotation. This polarized component due 
to E modes, which is ~ 3.5 yuK RMS at small angular scales 
plus an additional large-scale contribution of the North Galactic 
Spur of amplitude ~ 3 yuK, was then subtracted from the time- 
ordered data. Eurther iterations using the cleaned timelines were 
found to make a negligible difference. 

5.2. Results 

To make maps of the dipole, a second DaCapo run was made 
in the unconstrained mode for each LEI 70 GHz detector us- 


http://healpix.sourceforge.net 
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ing the polarization-cleaned timelines and the 30 GHz Madam 
mask, which allows 78% of the sky to be used. The extraction 
of the dipole parameters (Galactic latitude, longitude, and tem- 
perture amplitude) in the presence of foregrounds was achieved 
with a simple Markov chain Monte Carlo (MCMC) template¬ 
fitting scheme. Single detector hit maps, together with the white 
noise in the LFI reduced instrument model (RIMO), were used 
to create the variance maps needed to construct the likelihood 
estimator for the MCMC samples. Commander maps (Planck 
Collaboration X 2015) were used for synchrotron, free-free, and 
thermal dust for the template maps with the MCMC fitting for 
the best amplitude scaling factor to clean the dipole maps. The 
marginalized distribution of the sample chains between the 16th 
and 84th percentiles were used to estimate the statistical errors, 
which were 0?004, 0?009, and 0.16/iK for latitude, longitude, 
and amplitude respectively. The 50 % point was taken as the best 
parameter value, as shown in Table 1 . To estimate the system¬ 
atic errors on the amplitude in calibration process due to white 
noise, 1 // noise, gain fiuctuations, and ADC corrections, sim¬ 
ulated time-ordered data were generated with these systematic 
effects included. These simulated timelines were then calibrated 
by DaCapo in the same way as the data. The standard deviation 
of the input to output gains were taken as the error in absolute 
calibration with an average value of 0.11 %. 

Plots of the dipole amplitudes with these errors are shown 
in Fig. 6, together with the error ellipses for the dipole direc¬ 
tion. As can be seen, the scatter is greater than the statistical 
error. Therefore, we take a conservative limit by marginalizing 
over all the MCMC samples for all the detectors, which results 
in an error ellipse (+0?02, +0?05) centred on Galactic latitude 
and longitude (48?26, 264?01). The dipole amplitudes exhibit 
a trend in focal plane position, which is likely due to residual, 
unaccounted-for power in far sidelobes, which would be sym¬ 
metrical between horn pairs. These residuals, interacting with 
the Solar dipole, would behave like an orbital dipole, but in op¬ 
posite ways on either side of the focal plane. This residual there¬ 
fore cancels out to first order in each pair of symmetric horns in 
the focal plane, i.e., horns 18 with 23, 19 with 22, and 20 with 
21 (see inset in Fig. 6). As the overall dipole at 70 GHz is calcu¬ 
lated combining all the horns, the residual effect of far sidelobes 
is reduced. 


6. Validation of the caiibration and accuracy 
assessment 

In this section we present the results of a set of checks we have 
run on the data that comprise this new Planck release. Table 2 
quantifies the uncertainties that affect the calibration of the LFI 
radiometers. 

6.1. Absolute calibration 

In this section we provide an assessment of the change in the ab¬ 
solute level of the calibration since the first Planck data release, 
in terms of its impact on the maps and power spectra. Generally 
speaking, a change in the average value of G in Eq. (1) of the 
form 

<G) ^ <G) (1 + Sg), with Sg « 1 (40) 

leads to a change {T) (T) (1 -r ^g) in the average value of 

the pixel temperature T, and to a change ^ C^(l + 2Sg) in 
the average level of the measured power spectrum C^, before the 
application of any window function. Our aim is to quantify the 


Table 1. Dipole characterization from 70 GHz radiometers. 


Radiometer 

Amplitude 

Galactic Coordinates [deg] 

1 

b 

18M. 

. 3371.89 + 0.15 

264?014 + 0?008 

48?268 + 0?004 

18S. 

. 3373.03 + 0.15 

263?998 + 0?008 

48?260 + 0?004 

19M. 

. 3368.02 + 0.17 

263?981 +0?009 

48?262 + 0?004 

19S. 

. 3366.80 + 0.16 

264?019 + 0?009 

48?262 + 0?004 

20M. 

. 3374.08 + 0.17 

264?000 + 0?010 

48?264 + 0?005 

20S. 

. 3361.75 + 0.17 

263?979 + 0?010 

48?257 + 0?005 

21M. 

. 3366.96 + 0.16 

264?008 + 0?008 

48?262 + 0?004 

21S. 

. 3364.19 + 0.16 

264?022 + 0?009 

48?266 + 0?004 

22M. 

. 3366.61 +0.14 

264?014 + 0?008 

48?266 + 0?004 

22S. 

. 3362.09 + 0.16 

264?013 + 0?009 

48?264 + 0?004 

23M. 

. 3354.17 + 0.16 

264?027 + 0?009 

48?266 + 0?004 

23S. 

. 3358.55 + 0.18 

263?989 + 0?009 

48?268 + 0?004 

Statistical . . 

. 3365.87 + 0.05 

264?006 + 0?003 

48?264 + 0?001 

Systematic. . 

. 3365.5 + 3.0 

264?01 + 0?05 

48?26 + 0?02 

NominaP. . . 

. 3364.5 + 2.0 

264?00 + 0?03 

48?24 + 0?02 


This estimate was produced combining the LFI and HFI dipoles, and 
it is the one used to calibrate the LFI data delivered in the 2015 data 
release. 


value of the variation Sg from the previous Planck-L¥l data re¬ 
lease to the current one. We have done this by comparing the 
level of power spectra in the S = 100-250 multipole range con¬ 
sistently with Call3. 

There have been several improvements in the calibration 
pipeline that have led to a change in the value of (G): 


1. the peak-to-peak temperature difference of the reference 
dipole D used in Eq. (1) has changed by -rO.27 % (see 
Sect. 5), because we now use the Solar dipole parame¬ 
ters calculated from our own Planck measurements (Planck 
Collaboration I 2015); 

2. in the same equation, we no longer convolve the dipole D 
with a beam B that is a delta function, but instead use the full 
profile of the beam over the sphere (see Sect. 2); 

3. the beam normalization has changed, since in this data re¬ 
lease B is such that (Planck Collaboration IV 2015) 



5(6>,^)dQ < 1; 


(41) 


4. the old dipole fitting code has been replaced with a more 
robust algorithm, DaCapo (see Sect. 4). 


Table 4 lists the impact of such effects on the amplitude 
of fiuctuations in the temperature (T) of the 22 LEI radiome¬ 
ter maps. The numbers in this table have been computed by re¬ 
running the calibration pipeline on all the 22 LEI radiometer data 
with the following setup: 


1. a pencil-beam approximation for B in Eq. (1) has been used, 
instead of the full An convolution (“Beam convolution” col¬ 
umn), with the impact of this change quantified by Eq. (16), 
and the comparison between the values predicted by this 
equation with the measured change in the a^m harmonic co¬ 
efficients shown in Eig. 8 (the agreement is excellent, better 
than 0.03 %); 

2. the old calibration code has been used instead of the DaCapo 
algorithm described in Sect. 4 (“Pipeline upgrades” column); 
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Table 2. Accuracy in the calibration of LFI data. 


Type of uncertainty 

30 GHz 

44 GHz 

70 GHz 

Solar dipole. 

_ 0.10% 

0.10% 

0.10% 

Spread among independent radiometers^. 

_ 0.25 % 

0.16% 

0.10% 

Overall error. 

_ 0.35 % 

0.26 % 

0.20% 


^ This is the discrepancy in the measurement of the height of the first peak in the TT spectrum (100 <i< 250), as described in Sect. 6.4. 


o 

00 



Galactic longitude I [° ] 


o 



Fig. 6. Dipole amplitudes and directions for different radiome¬ 
ters. Top: errors in the estimation of the Solar dipole direction are 
represented as ellipses. Bottom: estimates of the amplitude of the 
Solar dipole signal; the errors here are dominated by gain uncer¬ 
tainties. Inset: the linear trend (recalling that the numbering of 
horns is approximately from left to right in the focal plane with 
respect to the scan direction), most likely caused by a slight sym¬ 
metric sidelobe residual, is removed when we pair the 70 GHz 
horns. 
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Fig. 7. Top: estimate of the value of cpu (Eq. 6) for each LFI 
radiometer during the whole mission. The plot shows the median 
value of (pD over all the samples and the 25th and 75th percentiles 
(upper and lower bar). Such bars provide an idea of the range 
of variability of the quantity during the mission; they are not an 
error estimate of the quantity itself. Bottom: estimate of the value 
of 0sky (Eq* 8). The points and bars have the same meaning as 
in the plot above. Because of the smallness of the value for the 
44 and 70 GHz channels, the inset shows a zoom of their median 
values. The large bars for the 30 GHz channels are motivated by 
the coupling between the stronger foregrounds and the relatively 
large power falling in the sidelobes. 


3. the 5 * Loai term has not been removed, as in the discussion 
surrounding Eq. (11) (“Galactic sidelobe removal” column); 

4. the signal D used in Eq. (1) has been modelled using the 
dipole parameters published in Hinshaw et al. (2009), as was 
done in Cal 13 (“Reference dipole” column). 

We have measured the actual change in the absolute calibration 
level by considering the radiometric maps (i.e., maps produced 
using data from one radiometer) of this data release (indicated 
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Table 3. Optical parameters^ of the 22 LFI beams. 


Radiometer 

/si [%] 

/d [%] 

0sky [^] 



70 GHz 



18M. 

0.38 

yj.yjy/ 

0.0000 

+0.0057 

-0.0058 

18S. 

0.62 

u.iUU _0 016 

0.0001 

+0.0102 

-0.0099 

19M. 

0.60 

0 140 +0003 
-0.007 

0.0001 

+0.0078 

-0.0077 

19S. 

0.58 

0 167 

-0.010 

0.0001 

+0.0083 

-0.0080 

20M. 

0.63 

0 157 

-0.003 

0.0001 

+0.0082 

-0.0079 

20S. 

0.70 

0 104 +0003 
-0.006 

-0.0000 

+0.0084 

-0.0085 

21M. 

0.59 

0 153 

_o.ooi 

-0.0000 

+0.0098 

-0.0100 

21S. 

0.70 

0 104 +0.007 
-0.003 

0.0000 

+0.0093 

-0.0090 

22M. 

0.44 

0 130 

_0 002 

0.0000 

+0.0074 

-0.0073 

22S. 

0.50 

0 167 

-0.004 

0.0000 

+0.0076 

-0.0078 

23M. 

0.35 

0 092 

yj.yjyz. _QQQ4 

0.0000 

+0.0054 

-0.0056 

23S. 

0.43 

0 122 

U.iZZ Q QQg 

-0.0000 

+0.0062 

-0.0061 



44 GHz 



24M. 

0.15 


0.0000 

+0.0059 

-0.0058 

24S. 

0.15 

/ _o.oooo 

0.0002 

+0.0073 

-0.0073 

25M. 

0.08 

0 0262+^-^^^^ 
^•^^'^^-0.0004 

0.0001 

+0.0040 

-0.0040 

25S. 

0.06 

o 

b 

o^ 

o o 

88 
o o 

0.0001 

+0.0032 

-0.0032 

26M. 

0.08 

88 
88 
d d 
+ 1 

(N 

O 

b 

0.0001 

+0.0038 

-0.0037 

26S. 

0.05 

o 

b 

o 

o o 

88 
o o 

0.0000 

+0.0030 

-0.0030 



30 GHz 



27M. 

0.64 

0 155 

_0.006 

0.0090 

+0.1475 

-0.1439 

27S. 

0.76 

0 190 

yj.Lyyj _ooo7 

0.0098 

+0.1656 

-0.1644 

28M. 

0.62 

0 1 S4 +0.005 
-0.012 

0.0063 

+0.1325 

-0.1310 

28S. 

0.83 

0 190 

yj.iyyj _ooi9 

0.0123 

+0.1810 

-0.1762 


^ The values for (pD and 0sky are the medians computed over the whole 
mission. Upper and lower bounds provide the distance from the 25th 
and 75th percentiles and are meant to estimate the range of variability 
of the quantity over the whole mission; they are not to be interpreted 
as error bars. We do not provide estimates for as they can all be 
considered equal to zero. 

with a prime) and of the previous data release and averaging the 
ratio: 

if = {^) - 1. m 

\ ^xxy I y 

where y indicates a radiometer at the same frequency as x and 
x', such that y x. The average is meant to be taken over all 
the possible choices for y (thus 11 choices for 70 GHz radiome¬ 
ters, 5 for 44 GHz, and 3 for 30 GHz) in the multipole range 
100 < T’ < 250. The way that cross-spectra are used in Eq. (42) 
ensures that the result does not depend on the white noise level. 
Of course, this result quantifies the ratio between the tempera¬ 
ture fiuctuations (more correctly, between the aim coefficients of 
the expansion of the temperature map in spherical harmonics) in 
the two data releases, and not between the power spectra. 

The column labeled “Estimated change” in Table 4 contains 
a simple combination of all the numbers in the table: 

Estimated change = (1 + 6beam) (1 + ^pipeline) x 

(l+6Gal)(l+6z))-l, ^ ^ 


^C/)^C/D^C/)^C/D^C/)^C/D^C/)^C/)^C/D^C/)^C/D 

(X)C»a)a)00-.--.-CMCMCOCO'^'^lOlOCDCDI^I^C»C» 

-.-■.-■.-■.-CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCM 



Fig. 8. Top\ comparison between the measured and estimated 
ratios of the aim harmonic coefficients for the nominal maps 
(produced using the full knowledge of the beam B over the 
An sphere) and the maps produced under the assumption of a 
pencil beam. The estimate has been computed using Eq. (11). 
Bottom: difference between the measured ratio and the estimate. 
The agreement is better than 0.03 % for all 22 LEI radiometers. 


where the e factors are the numbers shown in the same table and 
discussed above. This formula assumes that all the effects are 
mutually independent. This is of course an approximation; how¬ 
ever, the comparison between this estimate and the measured 
value (obtained by applying Eq. 42 to the 2013 and 2015 release 
maps) can give an idea of the amount of interplay of such effects 
in producing the observed shift in temperature. Eig. 9 shows a 
plot of the contributions discussed above, as well as a visual 
comparison between the measured change in the temperature 
and the estimate from Eq. (43). 


6.2. Noise in dipoie fitting 

We have performed a number of simulations that quantify the 
impact of white noise in the data on the estimation of the cali¬ 
bration constant, as well as the ability of our calibration code to 
retrieve the true value of the calibration constants. Details of this 
analysis are described in Planck Collaboration III (2015). We 
do not include such errors as an additional element in Table 2, as 
the statistical error is already included in the row “Spread among 
independent radiometers.” 
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Table 4. Changes in the calibration level between this (2015) Planck-LFl data release and the previous (2013) one. 



Beam 

Pipeline 

Galactic 

Reference 

Estimated 

Measured 

Radiometer 

convolution [%] 

upgrades [%] 

sidelobe removal [%] 

dipole [%] 

change [%] 

change[%] 




70 GHz 




18M. 

-0.277 

0.16 

0.000 

0.271 

0.15 

0.19 

18S. 

-0.438 

0.21 

0.000 

0.271 

0.04 

0.28 

19M. 

-0.443 

0.16 

0.000 

0.271 

-0.01 

0.09 

19S. 

-0.398 

0.21 

-0.002 

0.271 

0.08 

0.20 

20M. 

-0.469 

0.18 

0.000 

0.271 

-0.02 

0.18 

20S. 

-0.506 

0.19 

0.000 

0.271 

-0.05 

0.09 

21M. 

-0.439 

0.14 

0.000 

0.271 

-0.03 

0.01 

21S. 

-0.510 

0.21 

0.000 

0.271 

-0.03 

-0.30 

22M. 

-0.328 

0.20 

0.000 

0.271 

0.14 

0.08 

22S. 

-0.352 

0.20 

0.000 

0.271 

0.12 

0.21 

23M. 

-0.269 

0.32 

0.000 

0.271 

0.32 

-0.03 

23S. 

-0.320 

0.35 

0.000 

0.271 

0.30 

0.44 




44 GHz 




24M. 

-0.110 

0.28 

0.000 

0.271 

0.44 

0.71 

24S. 

-0.101 

0.26 

0.000 

0.271 

0.43 

0.34 

25M. 

-0.046 

0.30 

0.000 

0.271 

0.72 

0.27 

25S. 

-0.055 

0.31 

0.000 

0.271 

0.74 

0.62 

26M. 

-0.054 

0.19 

0.000 

0.271 

0.63 

0.52 

26S. 

-0.033 

0.17 

0.000 

0.271 

0.56 

0.44 




30 GHz 




27M. 

-0.498 

0.24 

-0.015 

0.271 

-0.01 

-0.15 

27S. 

-0.583 

0.14 

-0.019 

0.271 

-0.19 

-0.56 

28M. 

-0.440 

0.32 

-0.003 

0.271 

0.15 

0.35 

28S. 

-0.601 

0.32 

-0.004 

0.271 

-0.02 

0.22 

uncertainties 



In this data release, we have switched to the fr 


As discussed in Planck Collaboration IV (2015), the beams B 
used in the LFI pipeline are very similar to those presented 
in Planck Collaboration IV (2014); they are computed with 
GRASP, properly smeared to take into account the satellite mo¬ 
tion. Simulations have been performed using the optical model 
described in Planck Collaboration IV (2014), which was derived 
from the Planck radio frequency Flight Model (Tauber et al. 
2010) by varying some optical parameters within the nominal 
tolerances expected from the thermoelastic model, in order to 
reproduce the measurements of the LFI main beams from seven 
Jupiter transits. This is the same procedure adopted in the 2013 
release (Planck Collaboration IV 2014); however, unlike the case 
presented in Planck Collaboration IV (2014), a different beam 
normalization is introduced here to properly take into account 
the actual power entering the main beam (typically about 99 % 
of the total power). This is discussed in more detail in Planck 
Collaboration IV (2015). 

Given the broad usage of beam shapes B in the current LFI 
calibration pipeline, it is extremely important to assess their ac¬ 
curacy and how errors in B propagate down to the estimate of 
the calibration constants K in Eq. (1). 

In the previous data release we did not use our knowledge 
of the bandpasses of each radiometer to produce an in-band 
model of the beam shape, but instead estimated B by means of 
a monochromatic approximation (see Cal 13). In that case, we 
estimated the error induced in the calibration as the variation of 
the dipole signal when using either a monochromatic or a band- 
integrated beam, as we believe the latter to be a more realistic 
model. 


integrated beams produced using GRASP, which represents our 
best knowledge of the beam (Planck Collaboration IV 2015). 
We have tested the ability of DaCapo to retrieve the correct cal¬ 
ibration constants K for LFI19M (a 70 GHz radiometer) when 
the large-scale component (^ = 1) of the beam’s sidelobes is: (1) 
rotated arbitrarily by an angle -160° < 0 < 160°; or (2) scaled 
by +20 %. We have found that such variations alter the calibra¬ 
tion constants by approximately 0.1 %. However, we do not list 
such a number as an additional source of uncertainty in Table 2, 
since we believe that this is already captured by the scatter in the 
points shown in Fig. 10, which were used to produce the num¬ 
bers in the row “Inconsistencies among radiometers.” 

6.4. Inter-channel calibration consistency 

In this section we provide a quantitative estimate of the relative 
calibration error for the LFI frequency maps by measuring the 
consistency of the power spectra computed using data from one 
radiometer at time. By “relative error” we mean any error that 
is different among the radiometers, in contrast to an “absolute 
error,” which induces a common shift in the power spectrum. We 
have computed the power spectrum of single radiometer half¬ 
ring maps and have estimated the variation in the region around 
the first peak (100 < ^ < 250), since this is the multipole range 
with the best S/N. 

The result of this analysis is shown in Fig. 10, which plots 
the values of the quantity 


^rad 


(c™\ 


^ / rad 


(c?') 


1 , 


(44) 


freq 
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C»OOa)a)00-.--.-CMCMCOCO'^'^lOlOCDCDI^I^CX)CX5 

-.-■.-■.-■.-CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCM 



Fig. 9. Top: impact on the average value of the spherical har¬ 
monic coefficients (computed using Eq. 42, with 100 < ^ < 250) 
due to a number of improvements in the LFI calibration pipeline, 
from the first to the second data release. Bottom: measured 
change in the harmonic coefficients between the first and 
the second data release. No beam window function has been ap¬ 
plied. These values are compared with the estimates produced 
using Eq. (43), which assumes perfect independence among the 
effects. 



Fig. 10. Discrepancy among the radiometers of the same fre¬ 
quency in the height of the power spectrum Q near the first peak. 
For a discussion of how these values were computed, see the 
text. Inset: to better understand the linear trend in the 70 GHz 
radiometers, we have computed the weighted average between 
pairs of radiometers whose position in the focal plane is symmet¬ 
ric. The six points refer to the combinations 18M/23M, 18S/23S, 
19M/22M, 19S/22S, 20M/21M, and 20S/21S, respectively. Note 
that all six points are consistent with zero within Icr; see also 
Fig. 6. 


combinations of drad are consistent with zero within 1 cr (see the 
inset of Fig. 10). 

Since the cross-spectrum of two half-ring maps does not de¬ 
pend on the level of uncorrelated noise, the fluctuations of St 
around the average value that can be seen in Fig. 10 can be in¬ 
terpreted as relative calibration errors. If we limit our analysis to 
the multipole range 100 < ^ < 250, we can estimate the error of 
the 70 GHz map as the error on the average height of the peaks 
(i.e., the value cr/ ^^N, with cr being the standard deviation and 
N the number of points) that is, 0.25, 0.16, and 0.10 percent and 
30, 44, and 70 GHz, respectively. 


where is the cross-power spectrum computed using two 
half-ring maps, and (•) denotes an average over The quantity 
(Q)freq the Same average computed using the full frequency 
half-ring maps. Note that the dj-ad slope is symmetric around zero 
in the 70 GHz radiometers; this might be caused by residual un¬ 
accounted power in the far sidelobes of the beam. The same ex¬ 
planation was advanced in Sect. 5 to explain a similar effect. It 
is interesting to note that the amplitude of the two systematics 
is comparable; the trend in Fig. 6 has a peak-to-peak variation 
(in temperature) of about 0.5 %, while the trend in Fig. 10 has 
a variation (in power) of roughly 1.0 %. We combine the values 
of drad for those pairs of radiometers whose beam position in 
the focal plane is symmetric (e.g., 18M versus 23M, 18S versus 
23S, 19M versus 22M, etc.), since in these pairs the unaccounted 
power should be balanced. We have found that indeed all the six 


6.5. Inter-frequency calibration consistency 

In this section we carry on an analysis similar to the one pre¬ 
sented in Sect. 6.4, where we compare the absolute level of the 
maps at the three LFI frequencies, i.e., 30, 44, and 70 GHz. 

We make use of the full frequency maps, as well as the pair 
of half-ring maps at 70 GHz. Each half-ring map has been pro¬ 
duced using data from one of the two halves of each pointing 
period. We quantify the discrepancy between the 70 GHz map 
and another map by means of the quantity 


A 70 GHz,other 


^HRlxHR2 


^HR lx other 


- 1 , 


(45) 


where c'HR 1 xhr 2 cross-spectrum between the two 70 GHz 
half-ring maps, and cHRixother cross-spectrum between the 
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Fig. 11. Estimate of (Eq. 45), which quantifies the 

discrepancy between the level of the 70 GHz power spectrum 
and the level of another map. Top: comparison between the 
70 GHz map and the 30 GHz map in the range of multipoles 
100 < £ < 250. The error bars show the rms of the ratio within 
each bin of width 15. Bottom: the same comparison done be¬ 
tween the 70 GHz map and the 44 GHz map. A 60 % mask was 
applied before computing the spectra. 


first 70 GHz half-ring map and the map under analysis. In the 
ideal case (perfect correspondence between the spectrum of the 
70 GHz map and the other map) we expect = 0. As was the 
case for Eq. (44), this formula has the advantage of discarding 
the white noise level of the spectrum by using the cross¬ 
spectrum with the 70 GHz map, whose noise should be uncorre¬ 
lated. 

Over the multipole range 100 < 7” < 250, the average dis¬ 
crepancy^^ is 0.15 + 0.17% for the 44 GHz map, and 0.15 + 
0.26 % for the 30 GHz map, as shown in Eig. 11. Such numbers 
are consistent with the calibration errors provided in Sect. 6.4. 

6.6. Null tests 

In Cal 13 we provided a study of a number of null tests, with 
the purpose of testing the quality of the calibration. In this new 
data release, we have moved the bulk of the discussion to Planck 
Collaboration III (2015). We just show here one example, which 
is particularly relevant in the context of the LEI calibration vali¬ 
dation. Eigure 12 shows the variation in the quality of the maps 
due to the use of the full 471 convolution versus a pencil beam 
approximation, as discussed in Sect. 2. The analysis of many 
similar difference maps has provided us with sufficient evidence 
that using the full 471 beam convolution reduces the level of sys¬ 
tematic effects in the LEI maps. 


To reduce the impact of the Galactic signal we have masked 60 % of 
the sky, since we found that less aggressive masks produced significant 
biases in the ratios. 


7. Measuring the brightness temperature of Jupiter 

The analysis of the flux densities of planets for this Planck data 
release has been considerably extended. We only use Jupiter data 
for planet calibration, so we will focus the discussion on obser¬ 
vations of this planet. The new analysis includes all seven tran¬ 
sits of Jupiter through each main beam of the 22 LEI radiome¬ 
ters. The analysis pipeline has been improved considerably by 
taking into account several effects not included in Call3. 

Planets provide a useful calibration cross-check; in partic¬ 
ular, the measurement of the brightness temperature of Jupiter 
can be a good way to assess the accuracy of the calibration, as 
Jupiter is a remarkably bright source with a S/N per scan as high 
as 50 and a relatively well known spectrum. Eurthermore, at the 
resolution of LEI beams it can be considered a point-like source. 

7.1. Input data 

Table 5 lists the epochs when the LEI main beams crossed 
Jupiter, and Eigs. 13 and 14 give a visual timeline of these events. 
The first four transits occurred in nominal scan mode (spin shift 
2', 1° per day) with a phase angle of 340°, and the last three scans 
in deep mode (shift of the spin axis between rings of 0.5', 15' 
per day) with a phase angle of 250° (see Planck Collaboration I 
2015). The analysis follows the procedure outlined in Call3, but 
with a number of improvements: 

1. the brightness of Jupiter has been extracted from timelines 
to fully exploit the time dependence in the data; 

2. seven transits have been considered instead of two, which 
allowed us to better analyse the sources of scatter among the 
measurements; 

3. all the data have been calibrated simultaneously; 

4. different extraction methods have been exploited, in order to 
find the most reliable among them. 

In the following discussion, we refer to a “timeline” (one 
for each of the 22 LEI radiometers) as the list of values 
(7JCp,^, ATant,?), with t the epoch of observation, the 
instantaneous apparent planet positions as seen from Planck, Xt 
and if/t the corresponding beam pointing directions and orienta¬ 
tions, and ATant,? is the measured antenna temperature. The val¬ 
ues of ATant,? provided by the LEI pipeline are calibrated and 
have their dipole and quadrupole signals removed. The pipeline 
also provides the values of Xt and ij/t. We recovered JCp,^ from the 
Horizonson-line service. 

Samples from each radiometer timeline have been used in 
this analysis only if the following conditions were met: (1) the 
samples have been acquired in stable conditions during a point¬ 
ing period (Planck Collaboration II 2015); (2) the pipeline has 
not flagged them as “bad”; (3) their angular distance from the 
planet position at the time of the measurement is less than 5°; 
and (4) they are not affected by any anomaly or relevant back¬ 
ground source. We checked the last condition by visually in¬ 
specting small coadded maps of the selected samples. 

7.2. Description of the analysis pipeline 

In the following paragraphs, we describe how we improved the 
pipeline used to extract the brightness temperature Tb of Jupiter 
from the raw LEI data. Such extraction goes through a first esti¬ 
mation of the antenna temperature Ta and a number of correc¬ 
tions to take into account various systematic effects. We present 

http://ssd.jpl.nasa.gov/?horizons 
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Fig. 12. Difference in the application of the full An beam model versus the pencil beam approximation. Panel A \ difference between 
survey 1 and survey 2 for a 30 GHz radiometer (LFI-27S) with the An model, smoothed to 15°. We do not show the same difference 
with pencil beam approximation, as it would appear indistinguishable from the An map. Panel B\ double difference between the An 
1-2 survey difference map in panel A and the pencil difference map (not shown here). This map shows what changes when one 
drops the pencil approximation and uses the full shape of the beam in the calibration. Panel C: zoom on the blue spot visible at the 
top of the map in panel A. Panel D\ same zoom for the pencil approximation map. The comparison between panel C and D shows 
that the An calibration produces better results. 

Table 5. Visibility epochs of Jupiter. 


30 GHz 

44 GHz" 

70 GHz 

31 October-2 November 2009 

24-27 October 2009 

29 October-1 November 2009 

30 June-3 July 2010 

31 October-2 November 2009 

1-5 July 2010 

14-18 December 2010 

30 June-2 July 2009 

12-16 December 2010 

1-4 August 2011 

8-12 July 2010 

2-10 August 2011 

31 August-7 September 2012 

5-8 December 2010 

5-11 September 2012 

21 February-1 March 2013 

15-18 December 2010 

I- 3 August 2011 

7-9 August 2011 

31 August-6 September 2012 

II- 16 September 2012 

7-12 February 2013 

23 February-1 March 2013 

15-24 February 2013 


^ The observation of Jupiter is more scattered in time for the 44 GHz radiometers because of their peculiar placement in the LFI focal plane. 


the methods used to estimate Ta in Sect. 7.2.1, and then in 
Sect. 7.2.2 we discuss the estimation of T^. Since the computa¬ 
tion of Tb requires an accurate estimate of the planet solid angle 
Dp, we discuss the computation of this factor in a dedicated part. 
Sect. 7.2.3. 


7.2.1. Estimation of the antenna temperature 


Following Call3 and Cremonese et al. (2002), the recovery of 
the instantaneous planet signal from a timeline is equivalent to 
the deconvolution of the planet shape from the beam pattern Bt 
at time t. Since the planet can be considered a point source, the 
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Days after launch 

91 270 456 636 807 993 1177 1358 1543 


SSI SS2 SS3 SS4 SS5 SS6 SS7 SS8 



30 GHz 
44 GHz 
70 GHz 


2009 Oct 01 2010 Oct 01 2011 Oct 01 2012 Oct 01 2013 Oct 01 


Fig. 13. Visual timeline of Jupiters’s crossings with LFI beams. Here “SS” lables sky surveys. 


Days after launch 



20090ct01 20100ct01 2011Oct01 20120ct01 20130ct01 


Fig. 14. Time dependence of the angle between Jupiter’s direction and the spin axis of the Planck spacecraft. The darker horizontal 
bar indicates the angular region of the 11 LFI beam axes, and the lighter bar is enlarged by +5°. 


most practical way is to assume 


suming monochromatic radiometers): 


ATant,? = Ta,p + b, (46) 

where Ta,p is the unknown planet antenna temperature, h the 
background, and the beam response for the planet at 

the time of observation. Of course, Bt depends on the relative 
position of the planet with respect to the beam, If a suitable 

beam model is available, Bt can be determined and can be 
recovered from least squares minimization. We use an elliptic 
Gaussian centred on the instantaneous pointing direction as a 
model for the beam, because it shows a very good match with the 
main beam of the GRASP model (Planck Collaboration IV 2015), 
with peak-to-peak discrepancies of the order of a few tenths of a 
percent (the importance of far sidelobes is negligible for a source 
as strong as Jupiter). To compute Bt, the pointings are rotated 
into the beam reference frame, since this allows for better control 
of the beam pattern reconstruction.^^ 


7.2.2. Estimation of the brightness temperature 


Tb 


= B\ 


Planck 



Qb ^^Planck 

dT 



(47) 


where ^pianck is Planck’s blackbody function, Qb and Qp the 
beam and planet solid angles, and Tcmb = 2.7255 K is the tem¬ 
perature of the CMB monopole. In this 2015 Planck data release, 
we have also introduced corrections to account for the bandpass. 

The accuracy of the determination is affected by confu¬ 
sion noise (i.e., noise caused by other structure in the maps), 
which we estimate from the standard deviation of samples taken 
between a radius of 1° and 1.5° (depending on the beam) and 
5° from the beam centre. These samples are masked for strong 
sources or other defects. Since background maps are not sub¬ 
tracted, the confusion noise is larger than the pure instrumental 
noise. However, since the histogram is well described by a nor¬ 
mal distribution, we have used error propagation^^ to assess the 
accuracy of Tb against the confusion noise. 

In the conversion of into through Eq. (46), the 
pipeline implements a number of small corrections: 


In Calls, we computed the brightness temperature from the 
antenna temperature La by means of the following formula (as- 


This is the opposite of Cal 13, which used the planet reference 
frame. 


We can quickly derive an order of magnitude for the size of the 
confusion noise effect in the estimation of Ta. Since the confusion noise 
is of the order of a few mK and the number N of samples within 2 
FWHM (the radius used for estimating Tb) is of the order of 10^-10"^, 
we can expect an accuracy of the order of 1 mK/ ^^N ^0.1 mK. 
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1. detector-to-detector differences in the beam solid angle Qb, 
accounting for +6 %, which is probably the most important 
effect; 

2. changes in the solid angle of the planet, Qp, due to the change 
of the ^u^iiox-Planck distance, which introduces a correction 
factor of up to 6.9 % percent; 

3. changes in the projected planet ellipticity, due to the planeto- 
centric latitude of the observer and the oblateness of the 
planet, to reduce observations as if they were made at 
Jupiter’s pole; 

4. blocking of background radiation by the planet, changing 
from about 0.7 % to 1.5 %, depending on the ratio Qp/Qb; 

5. a 0si correction, which accounts for the fraction of radiation 
not included in the main beam (about 0.2 %). 

7.2.3. Determination of the solid angle of Jupiter 

The solid angle Qp of Jupiter for a given planet-spacecraft dis¬ 
tance Ap and planeto-centric latitude dp is given by 

Qp((5p) = j gp(48) 

where Ap^i-ef is a fiducial planet-spacecraft distance (for Jupiter, 
Ap,ref ~ 5.2 AU), is the solid angle of the planet as seen 

from its pole at the fiducial distance, and ep (< 1) the ratio be¬ 
tween the polar and equatorial radii of the planet. 

Our new pipeline does not use the ellipticity, FWHM, and 
orientation of the elliptical beam parameters provided by Planck 
Collaboration IV (2015), as they were derived from a marginal¬ 
ization over Ta,p on the same Jupiter data used in this analysis. 
However, because of a degeneracy between Ta,p and the beam 
parameters, the results were sensitive to details of the fitting pro¬ 
cedure (up to about 1 %), such as the radius of the area being 
analysed and the minimization method. Therefore, we have first 
determined a new set of beam parameters for each transit and 
computed the weighted average (using a numerical minimiza¬ 
tion), and then we used such parameters to determine rA,p. 

7.3. Results 

We present here the results per radiometer and transit, and we 
discuss also how to combine the measurements performed using 
the 30, 44, and 70 GHz radiometers to obtain three estimates of 
Tb at these three nominal frequencies. 

7.3.1. Brightness temperatures per radiometer and transit 

In Call3, based on the first two transits, we noted that confusion 
noise alone was not able to account for variations in Tb found for 
different radiometers belonging to the same frequency channel. 
We therefore assumed the presence of some unidentified system- 
atics error, dominating the ultimate accuracy in the measures. 

The fact that residual systematics are more important than 
confusion noise and background is still true in the 2015 data re¬ 
lease. The values of Tb measured by radiometers in the same 
frequency channel have a spread of 0.6 %, 1.0 %, and 0.6 % of 
the average signal at 30, 44 and 70 GHz, respectively, and are 
not normally distributed (see Fig 15). These observed disper¬ 
sions are a factor of ~ 3 larger than the confusion noise, and 
cannot be ascribed to the background, whose effect is only 1 % 
of the observed scatter. The excess dispersion must be due to a 
small residual systematic effect such as pointing, beam model, 
or mismatch of center frequency. 





Fig. 15. Distribution of the values of Tb for Jupiter measured 
by the 22 LFI radiometers during each of the seven transits. The 
histogram has been produced using five bins per frequency. The 
lack of Gaussianity in the three distributions is evident. 

One possible cause of the observed dispersion in the bright¬ 
ness temperature is some systematic effect in the estimation of 
the beam parameters (see Sect. 7.2.3). 

Non-Gaussianities in the beam, as well as beam smearing 
were investigated by replacing the elliptical beam with band av¬ 
eraged beam maps derived from GRASP calculations. The results 
are consistent with the elliptical beam, with residuals of at most 
4 X 10-3 K in T’a.p- 

We computed Ta,p and again using an analytical ap¬ 
proximation based on the assumption of negligible background 
(which is quite a good approximation for Jupiter), and compared 
the results. The two methods agree at the level of 6 x 10“^ K for 
Ta,p, and at the level of 6 x 10“"^ K for Tb. 

7.3.2. Combination of the results and comparison with 
WMAP 

We now want to combine the measurements of the 22 LFI ra¬ 
diometers in order to have three estimates of Tb at the three LFI 
nominal frequencies, 30, 44, and 70 GHz. 

Such determination of Tb depends on proper knowledge of 
the central frequency, Vcen, fof each detector. This parameter is 
derived from the bandpasses, but these are not known exactly. It 
is possible to remove most of the differences among the 30 GHz 
and among the 44 GHz radiometers by slightly changing the Veen 
values of the radiometers by as little as +0.2 GHz. However, it is 
still not clear how the bandshapes would have to be modified to 
explain such changes in Veen- For this reason, we did not include 
such corrections for 30 GHz and 44 GHz data. 

At 70 GHz the situation is complicated even more by the fact 
that the Veen values of the radiometers are significantly spread 
over the channel bandwidth, so that each of the 70 GHz radiome¬ 
ters samples a slightly different portion of the Jupiter spectrum. 
Indeed, comparing at 70 GHz from one transit to another af¬ 
ter having ordered the radiometers for increasing Veen. it is pos¬ 
sible to see quite a significant correlation (see Fig. 16). Using a 
linear regression of the at 70 GHz against Veen, a “cleaned” 
list of Tb values was obtained. Their average is identical to the 
simple weighted average of Tb for all of the transits and of the 
detectors Tb = (171.558 ± 0.008) K. The standard deviation of 
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Fig. 16. Distribution of the values ofTij for Jupiter as a function of the central frequency Vcen of each radiometer. 


Table 6. Brightness temperature of Jupiter. 


r'cen 

Tb 


[GHz] 

[K] 


28.4. 

145.9 + 

0.9 

44.2. 

159.8 + 

1.4 

70.4. 

171.6 + 

1.0 


the whole set of samples reduces to just 2.7 %. Interestingly, the 
inferred slope dT^jdvctn = (0.2570 + 0.0058) K/GHz matches 
very well with the one from WMAP data, dT^ldvc^n = (0.243 + 
0.025) K/GHz, and the correlation between transits loses most 
of its statistical significance. We attempted the same test for the 
30 GHz and 44 GHz channels, but the spread in Vcen is too small 
to produce meaningful results. These results open the possibility 
of including the spectral slope directly as a free parameter of the 
fit in a future analysis. 

The result of this analysis produced the brightness tempera¬ 
tures Tb listed in Table 6. The errors are 0.6 %, 0.9 %, and 0.6 % 
at 30, 44, and 70 GHz, respectively. 

As in Call3, in Fig. 17 we compared Jupiter’s averaged 
in each LFI band with the spectrum provided by WMAP. The 
agreement is quite good, with a difference that does not exceed 
0.5%. In this comparison, we must note that WMAP and Planck- 
LFI are calibrated on slightly different dipoles, with Planck 
assuming an amplitude of 3364.5/iK (Planck Collaboration 
I 2015), while WMAP assumed an amplitude of 3355yuK 
(Hinshaw et al. 2009). In addition, WMAP central frequencies 
are different from those of Planck. So to make the comparison, 
the measures provided in Weiland et al. (2011) are scaled by 
1.00268 and linearly interpolated to the averaged Vcen of each 
LFI frequency channel. 



30 40 50 60 70 

Frequency [GHz] 


Fig. 17. Top: brightness temperature of Jupiter (Tbr) compared 
with the data from Weiland et al. (201 1), linearly rescaled in fre¬ 
quency to match LFFs central frequencies Vcen and corrected for 
difference between LFFs and WMAP’s dipoles. Bottom: devia¬ 
tion from unity of the ratio between LFFs estimate for Tbr and 
WMAP’s. The agreement is excellent among the three frequen¬ 
cies. 
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7.4. What does the analysis of Jupiter tells about the LFI 
calibration? 

The errors in our estimates of Jupiter’s brightness temperature 
are 0.6%, 0.9%, and 0.6%. These are consistent with (albeit 
slightly larger than) the numbers presented in Sects. 6.4 and 

6.5, where our analysis on the consistency of inter-channel and 
inter-frequency cross-spectra in the range 100 < L < 250 pro¬ 
duced estimates between 0.1 % and 0.3 %. We believe that the 
larger error bars are a direct consequence of a number of facts: 
(1) bandpass uncertainties play a larger role in the analysis of 
a non-thermodynamic sky signal like Jupiter’s emission in the 
microwave range; (2) the amount of data usable for this analysis 
is only a limited fraction of the overall amount used to create 
LFI maps; and (3) uncertainties in the main beam profile have a 
larger impact in the study of a point-like source such as Jupiter, 
with respect to the analysis presented in Sect. 6, which dealt with 
angular scales corresponding to 100 < ^ < 250. 

8. Conclusions 

We have described the method used to calibrate the Planck/LFl 
data for the 2015 Planck data release, and we have provided a 
quantitative analysis which shows in detail the amount of change 
in the calibration level due to every improvement we imple¬ 
mented in the pipeline since the previous data release. Compared 
to our 2013 release, we have improved the LFI calibration in 
several ways, most notably in the use of an internally consistent 
dipole signal as a calibrator, as well as more accurate data anal¬ 
ysis algorithms. As a result, we have improved by a factor of 3 
the accuracy of the LFI calibration, which for the 2015 release 
ranges from 0.20 % to 0.35 %, depending on the frequency. 

An important byproduct of our analysis is a novel estimate of 
the solar dipole signal; using LFI data only, we have estimated 
the amplitude of the dipole to be 3 365.6 + 3.0yuK and the direc¬ 
tion of its axis to be (/, 90° -b) = (264.01 ± 0.05°, 48.26 ± 0.02°) 
(Galactic coordinates). This result matches the numbers pro¬ 
vided by Hinshaw et al. (2009) within Icr. This slight difference 
in amplitude has the effect of a shift by ~ 0.3 % in the over¬ 
all level of the calibrated timelines. Together with the improved 
LFI calibration, we now find very good agreement between the 
level of the spectra estimated by WMAR 
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Appendix A: Smoothing caiibration curves 

In this appendix we explain in some detail the algorithm used 
to smooth the LFI calibration constant timelines. This task is 
performed by the smoothing filter discussed in Sect. 3 (see also 
Fig. 3). 

A.1. Purpose of the smoother 

Variations in the orientation of the Planck spacecraft during the 
mission caused its instruments to observe the calibration signal 
(the CMB solar dipole) with varying amplitude. This in turns 
induced variations in the accuracy of the reconstruction of the 
calibration constant K = G~^ (Eq. 1), as shown in Fig. 4. Such 
variations are mainly due to statistical errors and are in general 
unrelated to the true stability of the LFI detectors. Therefore, 
before applying such calibration constants to the raw data mea¬ 
sured by the LFI radiometers, we have applied a low-pass filter 
that smooths most of the high-frequency fluctuations. 

In Cal 13 we employed a simple smoothing Alter, which used 
wavelets to clean the stream of calibration constants of high- 
frequency fluctuations (note that no smoothing was applied to 
the 30 GHz radiometers, since the calibration algorithm used for 
them made this step unnecessary). 

In the 2015 Planck data release, we have improved our 
smoothing Alter in order to take into account sudden jumps in 
the calibration constant that are caused by a genuine change in 
the state of the radiometer; we call these “real jumps.” Such real 
jumps are not due to statistical effects and it is therefore incor¬ 
rect to include such jumps when smoothing the data. The actual 
smoothing algorithm used in the 2015 Planck-L¥l data release 
works as follows: 

1. determine if there are sudden variations in the calibration 
constants that might be of non-statistical origin and make 
a list of them; 

2. split the stream of calibration constants into sub-streams, us¬ 
ing as boundaries the jumps found in the previous step; 

3. apply a low-pass Alter to each sub-stream defined in the pre¬ 
vious step. 

In the following paragraphs we will provide more details about 
the implementation of these steps. 

A.2. Detecting jumps in the caiibration constants 

We discuss here how our data analysis code is able to determine 
the presence of sudden jumps in the calibration constant that 
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are due to some real change in the state of the radiometer. This 
task is not trivial, since we must look for sudden variations in a 
stream of numbers (the K values) that is dominated by statistical 
noise; moreover, the rms of such data changes with time, since 
it is correlated with the amplitude of the dipole signal ATdip (see 
Fig. 4). Within a region of high rms, it is therefore possible to 
mistake a sudden change in the value of K due to statistical fluc¬ 
tuations with an intrinsic change in the radiometer’s calibration. 
We therefore have developed a flgure of merit which allows us 
to disentangle these two families of jumps. 

In order to implement our flgure of merit, we have deflned 
a procedure to quantify the level of statistical fluctuations in the 
data. Such a procedure is similar to a smoothing Alter and we 
describe it with the aid of Fig. A.l, which shows its application 
to some real data (the values of K for radiometer LFI-27M, cal¬ 
culated during the first two sky surveys). As already described 
in Sect. 3, the variation in the statistical noise in K (panel A) 
is related to the strength of the dipole signal, i.e., the amplitude 
of ATdip (panel B). Therefore, in order to properly weight the 
importance of variations in the value of K, we apply a moving 
average to the stream of K values, where the amplitude of the 
window (panel C) depends on the value of ATdip. The result of 
the moving average (A^sm) is shown in panel D of Fig. A.l. For 
each window, the code computes the rms, ctk, of the N values. 
This quantity depends both on the statistical noise and on the 
presence of real jumps in the value of K. Therefore, we use the 
expression 

ATdip X ^ (A.l) 

as a flgure of merit for the determination of the presence of real 
jumps. This quantity is proportional to the rms of the moving 
average, but it is weighted by the amplitude of the calibration 
signal; the stronger the latter, the more likely that a real jump is 
present in the N samples. 

The threshold used with Eq. (A.l) is deflned in terms of per¬ 
centiles, specifically, we consider all the data that are greater 
than the ^-th percentile to mark the presence of a real jump. The 
value for n depends on the radiometer: for 30 GHz radiometers 
it is 99; for 44 GHz it is 99.9; and for 70 GHz it is 99.5. Such 
values have been determined by considering the quality of the 
null^^ test (see Sect. A.4). 

A.3. The smoothing algorithm 

Once the positions of the jumps have been determined, the code 
applies a smoothing algorithm to each subset of values of the 
original K (produced by Da Capo) between two consecutive 
jumps. The algorithm applies a low-pass Alter in the Fourier do¬ 
main that retains only 5 % of the lowest frequencies. After this 
step, in order to further reduce the noise, we apply a moving av¬ 
erage to the result, where each sample is weighted by the value 
of ATdip. 

A.4. Validation of the algorithm 

Since an incorrect identification of a real jump is likely to pro¬ 
duce stripes in maps, we have run a number of null tests in order 
to optimize the free parameters of the smoothing algorithm (e.g., 
the threshold used to detect jumps and the widths of the moving 

The null tests we used in this process are typically the difference 
between survey maps; in the case of an ideal, perfectly calibrated in¬ 
strument, such differences should yield a map where all the pixels are 
zero. 



Fig. A.l. Example of application of the algorithm for the de¬ 
tection of jumps in the stream of ^ = G~^ values (Eq. 1) for 
LFI-27M during Planck's first year. Panel A: the set of calibra¬ 
tion constants K computed by Da Capo for LFI27M. The v-axis 
here is cropped to the first two sky surveys (one year of data). 
Panel B: the amplitude of the dipole signal, dT^ip, as seen within 
each ring by the Planck spacecraft. The two horizontal lines 
mark the thresholds for regions where the signal is considered 
either “strong” or “weak.” Panel C: the size N of the window 
used to compute the moving average of K (there are roughly 
A = 30 values of K per day). In regions where the dipole sig¬ 
nal is weak or strong, the window width is 1200 or 400 samples, 
respectively; outside such regions, we use a linear interpolation 
between these two values. Panel D: the result of applying a mov¬ 
ing average with the variable window size (Panel C) to the series 
of data shown in Panel A. This is not the smoothed series used 
for calibration; only the rms of the moving average is used (see 
next panel). Panel E: the flgure of merit used to detect jumps is 
the product of the dipole amplitude and the rms of the moving 
average. The threshold used to detect jumps in this particular ex¬ 
ample (LFI-27M) is equal to the 99th percentile of such values 
(grey dashed line). In this case, two jumps have been found (days 
257 and 450). 
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Fig. A. 2. Comparison between the behaviour of the smoothing code in the nominal case and in the case where the jump near day 
257 is considered to be a statistical fluctuation. Panel A\ value of the calibration constant in the two cases. Note that fast variations 
have been removed from the data in order to make the plot clearer. If the jump is not considered to be real by the algorithm, the 
smoothing stage introduces an increasing slope in the values (thin grey line). Panel B\ map of the difference between the maps from 
Surveys 1 and 2 (the jump happened during the first survey), in ecliptic coordinates. Residual systematic effects produce stripes that 
are either aligned with the direction of the scan (i.e., perpendicular to the ecliptic plane) or with the Galactic plane. The black square 
highlights a region in the map that has been observed during the jump near day 257. Panel C: zoom into the region highlighted in 
panel B (the difference between maps from Surveys 1 and 2), when the jump near day 257 has been considered real (see Panel A, 
thick black line). Panel D: the same as the previous panel, but data have been calibrated assuming no real jump near day 257 (Panel 
A, thin grey line). Features are sharper in the latter case, and therefore we can conclude that the former calibration produces better 
results. 


average windows). We have calibrated the data using a number 
of combinations of parameters and have produced single-survey 
sky maps (i.e., maps obtained using six months of data). We have 
then differenced them, under the hypothesis that a perfect cali¬ 
bration would produce a map where the value of every pixel is 
consistent with zero. Fig. A.2 shows an example of this analysis; 
this shows that the jump found by the code near day 257 is likely 
to be a real jump, because not considering it as such leads to a 
stronger stripe in survey-difference maps. 













